OIQP-06-14 



Effect of interaction shape on the condensed DNA toroid 
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PACS 82 . 35 . Lr - Physical properties of polymers 

PACS 87. lO.+e - General theory and mathematical aspects 

PACS 87 . 14 . Gg - DNA, RNA 

Abstract. - We investigate how different microscopic interactions between semiflexible chain seg- 
ments can qualitatively alter the physical properties of the condensed toroid. We propose a general 
form of the Hamiltonian of the toroid and discuss its analytic properties. For different interactions, 
the theory predicts different scaling behaviours of the mean toroidal and cross sectional radii, r c 
and r crosa , as functions of the contour length L: (r c , r cross) ~ L"^ Nc ^ with v = (|, |) for the van 
der Waals type, v = (—5, §) for the Coulomb type, v — ( — 1, 1) for the delta function type at- 
tractions in the asymptotic limit. For the toroids with finite winding number iV c = 100 ~ 400, we 
find v ~ for the Yukawa interaction with screening parameter k = 0.5 ~ 1.0, and v = 0.1 ~ 0.13 
for the van der Waals type interactions. These findings could provide possible explanation for 
the experimentally well known observation v ~ of the condensed DNA toroids. Conformational 
transitions are also discussed. 
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Condensed DNA toroid and its conformational proper- 
ties have been of great interest, for instance, as a possi- 
ble candidate for gene delivery in gene therapy [1,2]. Al- 
though condensed DNA toroids have been heavily investi- 
gated in experiment and theory ( [3-5] and the references 
therein), it remains unclear what physical factors deter- 
mine the toroidal sizes [2]. In fact, there is experimentally 
well known observation that the radius of DNA toroid r c 
is independent of chain contour length L = 400 — 50, 0006p 
(132.8 nm — 16.6 /im): r c ~ L v with v ~ [6]. In most an- 
alytic works, phenomenological free energy models balance 
the bending and surface free energies to estimate toroidal 
properties. The predicted scaling relation for the mean 
toroid radius is r c ~ L v with v = \ for large L [4,7,8]. 
This is however inconsistent to the experimental observa- 
tion. 

In this letter, we show how different microscopic inter- 
actions qualitatively modifies physical properties of the 
toroid, in particular, the mean toroidal and cross sectional 
radii. To achieve this, we first propose a general form 
of the Hamiltonian of the toroid and discuss its analytic 
properties. We then specify the explicit forms of attrac- 
tive interactions, namely delta-function, van der Waals, 
and Yukawa (screened Coulomb) type interactions. We 



show exponents of the mean toroidal and cross sectional 
radii in the large L (or equivalently large dominant wind- 
ing number N c ) limit are categorised into three distinct 
groups for different interactions. For experimentally re- 
alistic winding numbers N c = 100 ~ 400, we find v ~ 
for the Yukawa interaction with inverse screening length 
K = 0.5 ~ 1.0 and v = 0.1 ~ 0.13 for the van der Waals 
type interactions. We also describe conformational tran- 
sitions. 

Theory. — We propose the general Hamilto- 
nian of a tightly packed toroid, with any type of 
short/long-range attractive interactions between seg- 
ments: H ci (a, I, L, W) = fa 2 - W \?fV{N) + (L - 

Gap(N) . I denotes the persistence length of the semiflex- 
ible polymer chain of total contour length L. I is assumed 
to be large enough relative to the bond length lb to re- 
alise its stiffness (I ^S> lb)- a is an inverse toroidal radius 
and IF is a positive coupling constant of the two-body at- 
tractive interaction between polymer segments. N is the 
winding number: = [aL/27r]. 1 V(N) is the attractive 
potential in the toroid cross section in the unit of —W, 
which falls into the number of interacting segmental pairs 



( a ) present address. 



1 Gauss' symbol [x] gives the greatest integer among m < x. 
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for the van der Waals nearest neighbour interaction [5]: 
V(N) =3N- 2V3^/N-l/4. Gap(N) = V(N +1) - V(N) 
is introduced to compensate the continuity of the Hamil- 
tonian. Temperature j3 — 1/(&bT) is implicitly included 
in I and W. An important assumption is that the at- 
tractive potential (i.e. the second term in Ti. c i) is a 
linear function of the toroid radius, or equivalently of 
the mean perimeter of the toroid [9]. With the con- 
formation parameter c = ■frf'ji) > 0, and a variable 
x = fp > (i.e. [x] = N), the Hamiltonian is simplified 
by ttJ(o, I, L, W) = WL ■ H(c, x) and: 



H(c, x) 



/([*] 



4c 



Gap([x\), 



(1) 



where 



f(N) = NV(N + 1) - (N+l)V(N). 



For consistency, V(0) = V(l) = 0, so that -> at 

x = 0. By varying the value of W, V(2) can be normalised 
to unity. For a chain to form a toroid, VFL ^> 1 and c > 4 
are assumed. 

Most physical observables for the tightly packed toroids 
can be quite accurately estimated by the dominant toroid 
winding number N c . For instance, the mean toroid radius 
is given by 



2nN r 



L v(N a ) 



(2) 



with the exponent v(N c ). N c is given by the global min- 
imum of the Hamiltonian (1). Since TL(c,x) is governed 
by the value of c, our aim is now to derive the relation 
between c and iV c . 

Let us focus on the form of Tl(c, x) in the m-th segment 

(to < x < to + 1). The only extremum is given at x = 

i 

x c (m) = (2c/(to)) 3 . The condition for this existence in 
the segment is /(to) > and to < x c (m) < to + 1, from 
which it follows cl(to) < c < C[/(m), where cx(to) = 

2j£nj> cu{m) = 2f£l) ■ Provided f( m ) > 0, it means 
cu(m) ~ cl(to) for large to, and therefore, c ~ 2Wm) • 
Hence, for large N c , the c-N c relation is almost uniquely 
determined by 



2f(N c 



(3) 



Note that the second n-derivative of TL(c, x) is positive ex- 
cept at x S Z so that all existing extrema are stationary. It 
should be emphasised that the c-N c relation for a smaller 
value of N c is different from its asymptotic one. This be- 
comes crucial when we compare the theory to experiment 
in the end: real chains such as DNA have finite length. In- 
stead, we replace the above relation by c ~ CL ( N ^+ Cu ( Na "> ; 
which works well for the cases examined below. We now 
provide specific attractive potentials and calculate analyt- 
ical expressions for the mean toroidal and cross sectional 
radii. 




Fig. 1: A complete hexagon with a side of 5 segments in the 
toroidal cross section. 



The delta function potential: One of the simplest 
interaction is the delta functional attraction, ^^^(s) = 
—Wf^ds fuels' S(\r(s) — r(s')|) . In ref. [5], we have shown 
that the toroid of winding number N > 4 is the sta- 
ble ground state for c > 4. In this "ideal toroid" with 
zero thickness of the chain, every chain segment interacts 
equally with all the other segments accumulated on the 
same arc of the toroid. Hence, V(N) becomes the number 
of interacting segmental pairs: 

JV-l 

V(N) = N C 2 =J2 k = 2 N ^ N ~ 1) ' 

k=l 

Gap(N) is now given by Gap(N) = N. Consequently, 
f(N) is given by f(N) = ±N{N + 1) > 0. Substituting 
these into the expressions of cl,u{N), we obtain cl(N) = 

j^Xj- and cu(N) = ■ In the asymptotic limit, we 



have 



N c ~c, 



(4) 



N r and c 



as c 



dominant ideal toroid is 



u 



N c . Thus, the radius of the 



4tt/ 



2ttN c WL 



(5) 



Van der Waals potential: In real systems, the toroid 
has finite thickness hence physics could be different from 
that of the ideal toroid. Therefore, we must incorporate 
finite size effect into V(N) of the Hamiltonian, and then 
calculate the mean toroidal and cross sectional radii. 

First, we consider the van der Waals type interaction, 
or equivalently effective short-range dominant attraction. 
Toroidal cross section can be approximated by the hexag- 
onally arranged chains (Fig.l) [5], interacting via effective 
nearest neighbour van der Waals attraction. By definition, 
the excluded volume effect of the segments is incorporated 
in the hexagonally packed cross section. This hexagonal 
arrangement could be a good approximation as it has been 
experimentally observed for the condensed DNA toroid [2] . 

If the chains are packed in a complete hexagonal cross 
section (Fig.l), the winding numbers are N = 7,19,37, 
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Fig. 2: An exactly counted discrete function Vdiscrete(N) and 
its continuum approximation V(N) for the number of interact- 
ing pairs in the toroidal cross section. 



and so on. In such cases, the number of nearest neighbour 
interacting pairs between segments can be counted by the 
links between neighbouring pairs in the hexagonal cross 
section. We then obtain the number as a discrete function 
Vdiscrete(N). In the case of a complete hexagon with a 
side of (n + 1) segments, its winding number is given by 

n 

N= 1 + V6i = 3n(n + !) + !, 



i 

V?, 



N 



dominant toroid for large c: 



(8) 



Substituting this into r c = 2 ^ N , we obtain the mean 
toroid radius: 



(9) 



Note the mean toroid radius of T4 DNA in low ionic condi- 
tions and Sperm DNA packaged by protamines are quan- 
titatively fitted by this expression [5]. Also, the scaling 
property r c ~ Ls matches the one in [4,7,8]. 

Similarly, the mean radius of the toroidal cross section 
can be calculated for the complete hexagonal cross section 
with a side of (n + 1) monomers [5]: 

- 2W3/ 1\ 

Across — ^ 1^2/ 



2V3 + 3 A i 
12 



k, (10) 



When we connect the centres of nearest neighbour seg- 
ments in Fig.l, it results in a complete hexagon consisting 
of regular triangular cells with a side 1. Counting the 
number of the regular triangles' sides gives the number of 
the nearest neighbour interaction links: 

Vdi S crete{N) — [3 x (No. of the regular triangles : 6n 2 ) 

+ (Perimeter of the hexagon: 6n)] /2 
= (3 ■ 6n 2 + 6n) /2 = 3n(3n + 1) 
= 3N - 2\/3\/iV- 1/4, (6) 

for N = 1,7, 19, 37, ■ • •. Therefore, for a general value of 
N, we can approximate Vdiscrete(N) or analytically con- 
tinue it to the analytic function: 



where Id is the diameter of the segment. As N c ~ (2-^/3 c) 5 
for large c, we have 

3V3 + 6 r ||/wf 

r cross ^ ^2 ^ ' \~T ) 

Note that the scaling property r cross ~ hi is in agree- 
ment with the one in [4] obtained in the asymptotic 
limit. Also, we can formally consider the case of the 
ideal toroid (i.e. N c ~ c), although it has zero thickness: 

'cross — 2Atc 1j \21) td ' 



Yukawa potential and general theory. In exper- 
iment, when we put condensing agents such as multivalent 
cations into DNA solution, it can cause DNA to undergo 
the condensation from a worm-like chain (whip or coil) 
to toroidal states [1-6,8]. Due to surrounding ion clouds, 
effective interaction between the DNA segments could be 
described by the screened Coulomb (or Yukawa) potential: 



V(N) = 3N - 2V3^N-l/4. 



(7) 



V Y (r) 



-W 



e?q>(~Kl d (r/l d - 1)) 
(r/k) 



This approximates well an exactly counted discrete func- 
tion Vdiscrete(N) (see Fig. 2). Note that, up to JV = 3 we 
need not introduce the finite size effect, since there is no 
difference in the number of interacting pairs between the 
ideal and the van der Waals nearest neighbour toroids. 
Thus, N > 4 for this effect. Note also that the entan- 
glement (knotting) effect of the chain arrangement is ne- 
glected. 

The same analysis presented in the previous section 
leads to the following "asymptotic" c-N c relation of the 



where r is the distance between a pair of segments. The 
interaction range is characterised by the screening param- 
eter k (inverse screening length), which depends on salt 
concentrations. In the low screening limit k — > 0, the 
potential corresponds to the Coulomb interaction. 

Not only for the Yukawa potential with various n values, 
it is in general very difficult to analytically estimate the 
function V(N) just by counting the number of interacting 
segmental pairs. This is due to the long-range nature of 
the potential. Thus, we numerically compute the exact 
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value of the attractive potential in the toroidal cross sec- 
tion and do a fit by the following function. We assume 
that the function V(N) can be in general expanded as a 
polynomial in the radius of the cross section (oc (n + |)): 



V(.V) =lH,[n+^j +bi(n+ l - 



b 2 [n+- 



+b 3 [n 



-64 [n 



■ t)" J: ( a t) 2+ M n -\ 



+a 3 [n 



+A 4 



where A = ^0(^75) , 
A 3 = b 3 (^j: 3 A 4 = 



A 1 



A, 



(11) 



j_T" 



K 


a 


bo 







3.00007±0.00003 


8.38±0.00 


-6.25±0.02 


0.01 


2.65391±0.00529 


28.92±0.74 


-214.27±13.29 


0.1 


2.00845±0.00386 


110.45±2.02 


-954.14±37.64 


0.2 


2.00048±0.00207 


59.20±0.59 


-259.93±9.06 


0.3 


2.00008±0.00113 


41.06±0.22 


-121.61±3.17 


0.4 


2.00002±0.00063 


32.03±0.10 


-73.32±1.34 


0.5 


1.99957±0.00020 


26.69±0.02 


-51.12±0.25 


0.6 


1.99958±0.00013 


23.12±0.01 


-38.56±0.13 


0.7 


1.99963±0.00008 


20.58±0.01 


-30.68±0.08 


0.8 


1.99969±0.00006 


18.69±0.01 


-25.38±0.05 


0.9 


1.99974±0.00005 


17.23±0.00 


-21.63±0.04 


1.0 


1.99978±0.00004 


16.07±0.00 


-18.87±0.03 


13 = 6 


1.99999±0.00001 


9.56±0.00 


-6.85±0.00 



This fitting is based 
on the fact that V(N) of ideal and van der Waals types 
are given by V ldea i(N) = § (n + \) A - \) 2 - ^ = 

and Vvdw {N) = 9(n + \f - 6(n + \) + f = 
3iV — 2y/3^J N — j. All segments of the ideal toroid in- 
teract with each other with interaction energy — W, and 
V(2) is normalised to unity. This means the upper bound 
of V(N) should be Vi dea i(N). On the other hand, the lower 
bound of V(N) should be given by the van der Waals near- 
est neighbour toroid, interacting with energy — W. Hence, 
we have the inequality 



V ideal (N) > V(N) > V VDW (N). 



(12) 



Therefore, we speculate the Yukawa potential shall be in 
this region and might have a cubic term in (n + ^) in 
eq.(ll). For the analysis, the cross section has to be 
small relative to the mean toroid radius. Hence, a toroid- 
toroidal globule transition point (r c ~ r cross ) may get 
modified if we consider this effect more seriously. 
The corresponding function f(N) is given by: 

/(AO = NV(N +1) - (N + 1)V(Y) 



1 JVH-AfH 
16 



(13) 



where we sum over i = a, a— 1, a— 2, a— 3, a— 4 and C a = 
Aq, C a -i = A\, C a -2 = A2, C«_3 = A3, C Q _4 = A4. In 
the asymptotic limit but below the toroid-toroidal globule 
transition point, f(N) becomes Aq{% — ijN? for 2 < a < 

4, and A 1 (^--l)N 2 ^ for a = 2 (VDW type). Note that 
to satisfy the condition f(N) > (i.e. r c > 0), we require 
C a = A (~ 6 ) > for 2 < a < 4 and C a -i = Ax (~ 61) < 
for a = 2. 

The mean toroidal (r c ) and cross sectional (r cross ) radii 
are calculated for large c: 



Table 1: Least Squares Fitting of exactly computed V(JV) for 
the Yukawa potential with various screening parameter k and 
for — W( l -^Y with (3 = 6 using eq.(ll). For convenience, W 
and Id are taken to be unity. Data is presented for a, b and 
61 only. 

9 4ttZ 



26 WL ' ' c, ' oss 



I) Ideal type (a = 4): r c 
^( b f)Wl d L (6 >0), 

1 — - — 1 

II) Coulomb type (a = 3): r c = 67T^b 3 ( w ) 3 L~s , r croS s = 

^(^)%) § ^(^o>0), 
III) VDW type (a = 2): r c = (f , r cross = 

¥(^)Ww f (^<o). 

Note that the radii of the ideal and the van der Waals 
nearest neighbour toroids correspond to case I with bo = f 
and case III with b\ = —6, respectively. We will show 
below that case II is in fact the Coulomb type. 

Table 1 shows the least square fit of exactly computed 
V(N) for the potential — W(^f"f and the Yukawa poten- 
tial with inverse screening length k using eq.(ll). In 
numerics, we considered the complete hexagon with a 
side of up to 50 segments (i.e. winding number N = 
3n(n + 1) + 1 = 7351). The diameter of the chain l d 
is taken to be unity for brevity. We can confirm that 
the k = (Coulomb interaction) result is grouped by 
case II (Coulomb type). The Yukawa potential with 
k = 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 and (3 = 6 results are 
classified to case III (VDW type). Of great interest is the 
Yukawa potential with k = 0.01, 0.1, 0.2, where the expo- 
nent a takes the values 2^3. The transition from the 
Coulomb type (II) to the van der Waals type (III) occurs 
in this region. We note that for small n < 3, fit results 
deviate from V(N) of the Yukawa with n < 0.3. How- 
ever, it does not affect the large n behaviours, thus results 
reported here sustain. It should be also mentioned that 
the potential — W(^f) with (3 = 3 <~ 24 are categorised 
by case III (VDW type). But for /? = 1 ~ 2, cubic term 
becomes important (case II). 
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Fig. 3: The exponents u(N c ) of the mean toroid radius r c ~ 
l/( N c) f or ^ ne g n ^ e dominant winding number N c with screen- 
ing parameter k = 0.5 ~ 1.0, Coulomb, full and nearest neigh- 
bour Van der Waals, and delta function (ideal) attractions. 



The inequality V idea i{N) > V(N) > V V dw(N) means 
the exponents z^(A c ) for the mean toroid radius r c ~ 
[v(N c ) are b ounc ied by those of ideal and van der Waals 
toroids in the asymptotic limit: 



-l<u(N c ) < 



1 



(14) 



Also, we derive a significant fact that, as far as V(N) is 
given as a polynomial in (n + i) , we never have f(N c )^N c 
(i.e. N 2 ^c~L 2 ), thus r c ^L° for large c. This is con- 
tradicting to the experimentally well known observation 
r c ~ L . To resolve this problem, we plot the exponents 
v(N c ) of the mean toroid radius r c ~ L v ^ Nc ^ for the "fi- 
nite" dominant winding number iV c (Fig. 3). The expo- 
nents are now defined by u(N c ) = l — 2/i/f(N c ): 



r c = L/(2nN c ) = L/(2irc ) 



= (2tt) (W/21) "/ <Nc) L "/ <Nc) , 

where the finite function Vf(N c ) is defined as c = Nc f< " Nc \ 
We find that for N c = 100 ~ 400 (i.e. realistic winding 
number of DNA toroids such as T4 DNA [3] or Sperm 
DNA [6]), we have v ~ for the Yukawa interaction with 
k = 0.5 ~ 1.0, and v = 0.1 ~ 0.13 for the van der Waals 
interaction. These could explain the experimental obser- 
vation: v ~ 0. 

It is also possible to describe conformational transi- 
tions. When the radius of cross section becomes compara- 
ble to the segmental diameter r cross ~ l^, the whip-toroid 
(N c = 1) transition occurs. We have r c ~ L and the tran- 
sition line l/W ~ L 2 , which are independent of the inter- 
action shape. Moreover, for the toroid to toroidal globule 

2 i 

transition point r c ~ r cross , we have r c ~ l^L 3 . The scal- 
ing property r c ~ is similar to the expected scaling 
of globule like objects. Although the outer radius r outer 
scales differently for large c: r outer = r c + r cross — r cross , 
it scales the same as r c at the transition point. 



Conclusions. — To summarise, we have shown how 
different microscopic interactions between chain segments 
can alter the physical properties of the condensed toroid. 
In the asymptotic limit, exponents of the mean toroidal 
and cross sectional radii are categorised into three dis- 
tinct species: van der Waals type, Coulomb type, and ideal 
type. For the intermediate winding number of N c = 100 ~ 
400, we find v ~ for the Yukawa interaction with inverse 
screening length k = 0.5 ~ 1.0, and v = 0.1 ~ 0.13 for van 
der Waals interaction. These finding are consistent to the 
experimentally well known observation v ~ 0. It would 
be of great interest to check experimentally r c ~ L v for a 
fixed salt concentration, i.e., for a fixed screening parame- 
ter k. If we naively apply the asymptotic values of r c and 
f 'cross j transition lines are found to be interaction depen- 
dent and to agree with ref. [4] in VDW case. However, as 
has been shown, the asymptotic relations are not so accu- 
rate for finite and realistic winding number N c . Therefore, 
the transition lines are to be studied with a special care. 
Finally, It should be stressed that our generic theory can 
be applied straightforwardly to any toroidal condensation 
of semiflexible polymer chains with any type and number 
of microscopic interactions. 
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